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The release of the 2MASS Redshift Survey (2MRS), with it's 44599 galaxies 
allows the deduction of the galaxy's masses in nearly complete sample. A 
cubic box with side of 37 Mpc containing 2429 galaxies is extracted and the 
Newtonian gravitational field is evaluated both at the center of the box as 
well in 101 x 101 x 101 grid points of the box. The obtained results are then 
discussed at the hght of the shell theorem which states that at the internal 
of a sphere the gravitational field is zero, 
keywords 
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ture of the Universe 

1 Introduction 

The determination of the gravitational field in cosmology oscillates between 
the Newton law and various type of modifications to this law. The reference 
formula is the Newtonian force 



(1) 



where G is the gravitational constant, M the first mass, m the second mass 
and r the distance between the two masses. The enormous progresses in the 
observations of the spatial distribution of galaxies point toward a cellular 
structure, i.e. the galaxies are situated on the surfaces of bubbles, rather 
than to be aggregated in a random structures, see Coil (2012). In the limiting 
case in which all the galaxies were situated on the surfaces of spheres the 
gravitational forces should be zero due to the shell theorem or nearly zero 
due to the fact that the galaxies are distributed in a discrete way rather than 
in a continuous way. This paper describes in Section [2] two astronomical 
catalogs which allow to calibrate the size of the cosmic voids. Section [3] is 
devoted to the study of the photometric properties of of a nearly spherical 
distribution of galaxies and to a careful analysis of completeness connected 
with the selected astronomical catalog. Section |4] contains the evaluation of 
the Newtonian gravitational field in a box of 37 Mpc in which the boundary 
conditions are properly evaluated. Section [5] reports a comparison between 
three ideal structures and a real void as extracted from a slices oriented 
catalog. 



2 Observations 

This section processes the Sloan Digital Sky Survey Data Release 7 (SDSS 
DR7) 
see 



see 



Huchra et al. (2012) 



Abazajian et al. (2009), and the 2MASS Redshift Survey (2MRS), 



2.1 Observed statistics of the voids 

The distribution of the effective radius between the galaxies of SDSS DR7 



has been reported in Pan et al. (2011). This catalog contains 1054 voids and 



Table [T] reports their basic statistical parameters. 
2.2 The 2MASS 

The 2MASS is a catalog of galaxies which has instruments in the near-infrared 
J, H and K-bands (1-2.2 m) and therefore detects the galaxies in the so 
called "Zone of Avoidance," see iJarrettl (|2004l); iCrook et all (|2007[). At the 



moment of writing the 2MRS consists of 44599 galaxies with redshift in the 
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Table 1: The statistical parameters of the effective radius in SDSS DR7. 



parameter value 

mean 18.23/i~^ Mpc 

variance 23.32/i~^ Mpc^ 

standard deviation 4.82/i~-'^ Mpc 

kurtosis 0.038 

skewness 0.51 

maximum value 34.12/i~^ Mpc 

minimum value 9.9/i~^ Mpc 



Huchra et al. (2012). 



interval < ^ < 0.09, see 

galactic latitude, the galactic longitude and the expansion velocity 



The catalog contains the 
from 



these three parameters is possible to deduce the Cartesian coordinates, X, Y 
and Z expressed in Mpc. Figure [l] reports a cut of a given thickness of 2MRS 
where A express the thickness of the cut and No the number of selected 
galaxies. 



3 Photometric properties 



This section reviews the photometric maximum in the framework of the lumi- 
nosity function for galaxies and the Malmquist bias which fixes the concept 
of complete sample. A model for the luminosity of galaxies is the Schechter 
function, L*,a, $), where a sets the slope for low values of L, L* is the 
characteristic luminosity, and $* is a normalization, see (Zaninetti, 2010b 



Eqn. (55)). This function was suggested by Schechter (1976) and the distri 



bution in absolute magnitude ,$(M; M*, a, $), can be found in (Zaninetti 



2010b, Eqn. (56)) where M* is the characteristic magnitude as derived from 
the data. The parameters of the Schechter function concerning the 2MRS as 



well the bolometric luminosity, M|^^, can be found in Cole et al. (2001) and 
are reported in Table |2] The number of galaxies at a given f lux / as a func 
tion of the redshift z, ^ni^^j {z; Zcrit, c, Hq), can be found in (|Padmanabhan 



dfldzdf ^ 
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Figure 1: Cut of the 3D spatial distribution of 2MRS in the X = plane 
when A = 10 Mpc, the squared box has a side of 148 Mpc and Nq = 1244. 



Table 2: The parameters of the Schechter function and bolometric magni- 
tude for the 2MRS in the — band. 



parameter 2MRS 
M* - 5 logio h [mags] (-23.44 ± 0.03) 
a -0.96 ± 0.05 

[h^ Mpc-^] ((1.08 ± 0.16)10-2) 

3.39 

h 0.7 
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Figure 2: The galaxies of the 2MRS with 8.48 < m < 10.44 or 
1202409^^ < / < 7267112^ are organized in frequencies versus heho- 
centric redshift, (empty circles); the error bar is given by the square root of 
the frequency. The maximum frequency of observed galaxies is at z = 0.015. 
The full line is the theoretical curve generated by ^j^j^^j {z] Zcrit, c, Hq). In 
this plot M®^ = 3.39, h = 0.7, M* =-24.87, a =-0.98, $*=0.0037, = 721 
and the number of bins 40. 



1996 Eqn.(1.104)) or (Zaninetti 2010b, Eqn. (6)), where c/fi, dz, and df are 



the differentials of the solid angle, the red-shift, and the flux, respectively, 
Zcrit a parameter. Ho the Hubble constant and c is the velocity of light. The 
number of galaxies at a given flux has a maximum at z = Zmax{zcrit, see 



(Zaninetti, 2010b, Eqn. (8)). Figure ^reports the number of observed galax- 



ies in the 2MRS catalog at a given apparent magnitude and the theoretical 

dN , 
dfldzdf ^ 



curve as represented by {z; Zcrit, c, Hq). The merit function can be 



computed as 



= 1 '^riastriz) 



where n is number of data, the two indexes theo and astr stand for the- 
oretical and astronomical, respectively and (Jnastriz)'^ is the variance of the 
astronomical number of data; the obtained value is reported in the caption 
of Figure [2j 

The total number of galaxies in the 2MRS as function of z is reported 
in Figure [3] as well as the theoretical curve represented by the numerical 
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Figure 3: All the galaxies of the 2MRS with rriks < 11.75 are organized in 
frequencies versus heliocentric redshift, (empty circles); the error bar is given 
by the square root of the frequency. The maximum frequency of all observed 
galaxies is at z = 0.017. The full line is the theoretical curve generated by 
^^^{z; Zcrit, c, Ho). In this plot Mf^ = 3.39, h = 0.7, M*=-23.97, a =-0.96, 
$*=0.0037, = 1267 and the number of bins is 30. 



integration of 



dN 



dfldzdf * 

The mass of a galaxy can be evaluated once the mass luminosity ratio, R 
is given by 

R-if) . (3) 



Some values of R are now reported: i? < 20 by Kiang (1961) and Persic 



& Salucci 



der Marel 



(1992), i? = 20 by Padmanabhan (1996) and R 



(1991). The Malmquist bias, see Malmquist (1920 



originally applied to the stars and later on to the galaxies by Behr 



5 .93 b y 
1922|), 



van 



(1951). 



was 
The 



observable absolute magnitude, MLirriL] z, Hq), as a function of the limiting 
apparent magnitude, rriL, can be found in (Zaninetti, 2010b, Eqn. (51)). 



The bias predicts, from a theoretical point of view, an upper limit for the 
maximum absolute magnitude which can be observed in a catalog of galaxies 
characterized by a given limiting magnitude and Figure |4] reports such a 
curve as well the galaxies of the 2MRS. 

The limiting magnitude of the 2MRS is mL=11.75 and therefore the 
2MRS is complete for z < 0.0025. For values of z greater than this value 
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Figure 4: The absolute magnitude M of 36464 galaxies belonging to the 
2MRS when = 3.39 and Hq = 70 km s"^ Mpc~^ (green points). The 
upper theoretical curve as represented by ML^niL; z, Hq) is reported as the 
red thick line when 771^=11.75. 



the observed sample is not complete and we can introduce the efficiency, 
es{z;Mmax,Mmin,mL,c,h), where Mmax and Mmin are the maximum and 
minimum absolute magnitude of the considered catalog and h = ifg/lOO, 
see (Zaninetti, 2010b, Eqns. (51-53)). As an example when z ~ 0.017 the 
sample covers the 52% of the range in absolute magnitude. 



4 The gravitational field 

We now explore a connection with the shell theorem. Let us consider a 
spherically symmetric surface of radius a on which a total mass M is disposed 
in a uniform way. The mass per unit area /z is 



M 



The force at the internal of the spherical surface is 

r = r <a , 



(4) 



(5) 



see Eq. (11.24) in Alonso & Finn (1992) and therefore the shell theorem can 
be formulated "A uniform shell of matter exerts no gravitational force on a 
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particle situated inside a shell". We now compute the force on the center 
{x—0,y=0,z=0) of a hemisphere which resides on the positive z-axis. The 
vectorial intensity of the field is 

rfr = G'^ , (6) 

being 

dm = fia^dQ , (7) 

where dfl = smOdOdcj) is the solid angle with < < 27r and < 6* < |. 
The three Cartesian components of the field in the 3D case are 

d(f) G ji sinO cosO dO 

/■27r /-O 

Fx — J d(f) j G ji sin9 sinO cos d9 (8) 

Ty — d(f) G iJ, sinO sinO sin dd . 
Jo if 

The integration gives in the 3D case for the three forces at the center 

r. = r, = o . (9) 

We now consider the 2D case of mass concentrated on a half circle of radius 
a situated on the positive y-axis, where now the mass per unit length /X2D is 

M , , 

The 2D vectorial intensity of the field is 

^ = G% , (11) 

being 

dm — ii2D(id9 , (12) 
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where < < vr. The 2D Cartesian gravitational components of the force 
at the center {x=0,y=0) are 

r =G— r sindde 
a Jo 

= r cos Ode . (13) 

a Jo 

The integration of the 2D case gives 

^ 2GM 



r. = . (14) 

At the moment of writing the Committee on Data for Science and Technology 
(CODATA) recommends 

G = (6.67384 ± 0.00080) x 10""- — - , (15) 

kgs'^ 



see Mohr et al. (2008). Before to continue we express the Newtonian gravi- 
tational constant in the following units: length in Mpc, mass in Mgai which 
is 10^^M^,Q and yrS which are 10^ yr 

G = 4.4997510-^-^^^^ . (16) 



The two formulas (|9j) in 2D and (14) in 3D represent a useful reference to 
test a numerical code and to fix the range of variability of the gravitational 
field. According to our 3D theory the gravitational field at the center of the 
cosmic voids varies between the minimum value of zero (shell theorem) and 
a maximum value 

T,= ^^6.769W-^nMpE^ , (17) 

2R yr8^ 

where N is the number of galaxies in the spherical shell surrounding the 
cosmic void having mass M = Mgai and i?=18.23 Mpc is the average radius 
of the cosmic voids. 

We are now ready to process the 2MRS data and we associate to each 
galaxy, as reported in Figure[T]a mass given by Eqn. [3j The three components 
of the gravitational field are reported in table [3) 
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Table 3: 3D Gravitational forces expressed in at the center of a 3D 

box of side 37 x 2 Mpc when R =6 and theoretical 3D formula ( |9|. At 
z=0.008 the efficiency of the sample is ~ 70.6 %. 



Environment 








||r|| 


real structure 


-9.7710-6 


-1.5310-5 


-3.0410-6 


1.8410-5 


half sphere 








2.5510-5 


2.5510-5 



From a careful analysis of Table |3] it is possible to conclude that the 
gravitational field is greater than zero but smaller in respect to the case in 
which all the galaxies resides on a half sphere of radius equal to the averaged 
radius of the sample. Figure |5] reports a slice at the middle of a smaller box. 

The spatial displacement of the 3D grid r{i,j,k) which represent the 
absolute value of the gravitational field can be visualized through the iso- 
density contours, and as an example we considered a 101x101x101 grid. In 
order to do so, the maximum value ^n{i, j, k)max minimum value 

r„(i, j, should be extracted from the three-dimensional grid. A value 

of this grid can be fixed by the following equation: 

^nii y j y chosen ^nij' y j y k^min ~l~ (J^ni^i^ j j k^max ^ ni^i j j j k^min) ^ COCy , 

(18) 

where coef is a parameter comprised between and 1. This iso-surface render- 
ing of the gravitational field is reported in Fig |6j the Euler angles character- 
izing the point of view of the observer are also reported. Another interesting 
quantity to plot is the statistics of the values of the already defined spatial 
grid r„ which holds 101 x 101 xlOl values of gravitational field, see Figure 
[7j From this histogram is possible to conclude that the 90% of the space has 
a gravitational field comprised 3.24 10-^ ^^^'^g""' < r„ < 3.16 10'^^^^!^. 
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Mpc from center 



Figure 5: Cut-map of the 3D gravitational field of 2MRS when R=6. In 
order to have periodic boundary conditions the side of the box is 37 Mpc. 




Figure 6: Iso-surface of the logarithm of the 3D gravitational field of 2MRS 
when R=6 and coef = 0.43. The orientation of the figure is characterized 
by the Euler angles, which are $=30°, 6=30° and \1>=30°. 
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Figure 7: Decimal logarithmic histogram (step-diagram) of the values of the 
gravitational field evaluated in 101x101x101 points. 
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Figure 8: Polar plot of the real galaxies (green points) belonging to the sec- 
ond CFA2 redshift catalog and the simulated galaxies in the PVT framework 
(red points). More details can be found in Zaninetti (2006). 



5 The Voronoi simulation 

The Poisson Voronoi tessellation (PVT) is a useful tool to explore the spatial 
clustering of galaxies. The filaments of galaxies visible in the slices-type 
catalogs are due to the intersection between a plane and the PVT network of 
faces as a first approximation. An improvement can be obtained by coding 
the intersection between the slice of a given opening angle and the PVT 
network of faces, see Zaninetti (2006, 2010a). As an example Figure [s] reports 
both the CFA2 slice as well the simulated slice. 



We now test formula (14) in a discrete environment rather than in the 
continuous case. The test now calculates the two forces P^; and P^ in the 
center of the circle and in a 2D irregular Voronoi polygon generated by PVT 
which has the same averaged radius of the circle and center occurring in the 
same location of the generating seed. The half Voronoi polygon and the half 
circle are displayed in Figure |9] and the two forces P^; and Ty in Table [4j We 
are now ready to process a real void and our attention is focused on a CFA2 



slice visible in Figure 10 



A real void is extracted and the averaged radius of the galaxies on the 
boundary of that void is computed, see Figure [TT] The forces in the x and y 
direction are then computed and reported in Table |5} 

The presence of both a discrete number of galaxies and a non exactly 
symmetric displacement of the galaxies produces gravitational forces that 
take a finite value rather than zero. It is interesting to point out that Pj^ 
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Figure 9: Half circle with a=12.25 Mpc (green stars) and half irregular 
polygon (red squares). The number of unit masses, 1 Mgai is 455. 



Table 4: Gravitational forces expressed in 



yrS^ 



in the comparison be- 



tween half circle, 2D formula (14), and half irregular polygon. The parame- 
ters are Mgai=455 and a = 12.25453Mpc. 



Environment 






Half circle — theory 





8.6710-6 


Half circle — numeric 


-4.16 10-« 


8.81 10-6 


Half polygon — numeric 


1.3710-6 


1.3610-5 
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Figure 10: Polar plot of the real galaxies (green points) belonging to the 
second CFA2 redshift catalog. 



Table 5: Gravitational forces expressed in ^^ga''"' for the comparison be- 
tween a circle and a real CFA2 void. The parameters are Mgai=101 and 
a = 12.25453Mpc. 



Environment 






circle theory 








half circle theory 





1.9210-6 


real void numeric 


-1.36 10-^ 


8.310-^ 
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Figure 11: Circle with a= 12.25 Mpc (green stars) and real void extracted 
from a CFA2 slice (red squares). The number of galaxies with unit mass is 
101. 
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due to the galaxies on the boundary of the real void is smaller than the 
theoretical value as given by the half circle which represents a maximum 
theoretical value. 



6 Conclusions 

The masses of the galaxies can be deduced starting from the luminosities in 
the framework of the mass luminosity ratio, R. The spatial distribution of the 
masses of the galaxies allows the computation of the Newtonian gravitational 
forces on the unit mass. As a reference for the evaluation of the forces the 
2D and 3D shell theorem is analyzed. The evaluation of forces at the center 
of the box allows to conclude that the forces are smaller in respect to the 
mass concentrated on a half sphere of radius equal to the averaged radius 
of the selected sample of galaxies but bigger than zero due to the fact that 
the distribution of the galaxies is discrete rather than continuum. A careful 
analysis of a cubic box having side of 37 Mpc allows to say that the 90 % of the 
space has gravitational forces around the average value of 2.1 iQ-s -^^P^^g"' . 
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